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ABSTRACT 

In this work, we explore the effects of burst rate density evolution on the observed brightness distribution of 
cosmological gamma-ray bursts. Although the brightness distribution of gamma-ray bursts observed by the 
BATSE experiment has been shown to be consistent with a nonevolving source population observed to red- 
shifts of order unity, evolution of some form is likely to be present in the gamma-ray bursts. Additionally, 
nonevolving models place significant constraints on the range of observed burst luminosities, which are 
relaxed if evolution of the burst population is present. In this paper, three analytic forms of density evolution 
are examined. In general, forms of evolution with densities that increase monotonically with redshift require 
that the BATSE data correspond to bursts at larger redshifts, or to incorporate a wider range of burst lumi- 
nosities, or both. Independent estimates of the maximum observed redshift in the BATSE data and/or the 
range of luminosity from which a large fraction of the observed bursts are drawn therefore allow for con- 
straints to be placed on the amount of evolution that may be present in the burst population. Specifically, if 
recent measurements obtained from analysis of the BATSE duration distribution of the actual limiting redshift 
in the BATSE data at z lim = 2 are correct, the BATSE N(P) distribution in a A = 0 universe is inconsistent at 
a level of ~ 3 a with nonevolving gamma-ray bursts and some form of evolution in the population is required. 
The sense of this required source evolution is to provide a higher density, larger luminosities, or both with 
increasing redshift. 

Subject headings: cosmology: theory — gamma rays: bursts 


1. INTRODUCTION 

The paradigm of a cosmological origin for gamma-ray 
bursts has gained greater acceptance by the scientific com- 
munity in light of the Burst and Transient Source Experiment 
(BATSE) results which indicate that the gamma-ray bursts are 
distributed isotropically on the sky, yet possess an integral 
brightness distribution that deviates significantly from the 
— 3/2 power law indicative of homogeneity in Euclidean space 
(Meegan et al. 1992). Numerous models of cosmological burst 
production have been postulated, including NS-NS mergers 
e.g., Piran et al. (1991), emission from massive black holes in 
the center of galaxies (Brainerd 1994), and failed supernovae 
(Woosley 1993). Through a variety of analysis techniques, 
Wickramasinghe et al. (1993), Mao & Paczynski (1992), Piran 
(1992), Fenimore et al. (1993), and others have shown that the 
observed brightness distribution of gamma-ray bursts, 
jV(>P), is consistent with a population of mono-luminous 
sources distributed with a constant comoving rate density out 
to a limiting redshift of z max ~ 1 in a Friedmann universe with 
cosmological constant A = 0. 

Emslie & Horack (1994) have relaxed the restriction of 
mono-luminosity and also allow a general value for the cosmo- 
logical constant A in their analysis of the BATSE brightness 
distribution. Utilizing an integral moment analysis technique, 
they find not only that the limiting redshift for most cosmo- 
logical scenarios is indeed of order unity, but also that the 
gamma-ray burst luminosity function must be fairly narrow for 
such scenarios, with most observed bursts originating from a 
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luminosity range that is no wider than a factor of ~ 10. Emslie 
and Horack additionally show that if one assumes the 
observed luminosity function is substantially wider than this 
value, the observed BATSE brightness distribution is consis- 
tent with a constant comoving number density only in the 
context of an accelerating universe with a positive value for A. 

None of these previously mentioned analyses examine in 
detail the occurrence of evolution in either the comoving 
source density function n c or in the luminosity function 4>(L). 
Evolution, either of the form n c (z), <j>(L, z), or both may be 
present at some level. It is therefore important to explore the 
effects of evolution on the observed GRB brightness distribu- 
tion and to determine the consistency of various evolutionary 
scenarios with the observed N(P) distribution. 

In this work, we utilize the expanded dataset offered by the 
second BATSE gamma-ray burst catalog (Meegan et al. 1994). 
We continue to assume that the luminosity function <p{L) is 
independent of z, regardless of its particular form, but we do 
allow the comoving rate number density of gamma-ray bursts 
to be a function of redshift z [or some equivalent parameter 
such as cosmic look-back time t(z)]. This mode of evolution is 
known as pure density evolution (PDE), and has been adopted 
for several reasons. First, because we are only beginning to 
explore the effects of source evolution in the context of cosmo- 
logical models, it makes sense to begin with the simplest set of 
assumptions and grow into more complex evolutionary sce- 
narios only if necessary. PDE is the simplest mathematical 
representation of evolution, where all evolutionary effects are 
contained in the comoving density function n c (z). Second, 
because the luminosity function is independent of redshift, the 
effects of n c (z) and (j>(L) are separable and hence the method- 
ology of Emslie & Horack (1994), based on the technique of 
integral moment analyses, can be applied to the observed 
brightness distribution to determine the range of cosmological 
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scenarios that are compatible with the data. [Redshift-depen- 
dent luminosity functions that are separable, i.e., 
4>(L, z) = <p(L)g(z), can also be analyzed under this method- 
ology.] Third, PDE may indeed accurately describe the evolu- 
tionary situation for a number of cosmological burst models, 
for example the NS-NS merger scenario where the burst’s 
luminosity depends only on z-independent physical quantities 
such as the mass of a neutron star, but the event rate would 
depend on redshift because the formation of progenitor 
systems traces galactic star formation which is likely to vary in 
time (e.g., Hartmann 1995). 

We examine the compatibility between three different evolu- 
tionary scenarios and the brightness distribution of the second 
BATSE catalog of data (Meegan et al. 1994): 

1. Power-law evolution, where the comoving rate density of 
bursts is proportional to (1 + zf, where /? is a parameter. 

2. Exponential evolution, where n c (z ) is proportional to exp 
(rjxH 0 ), where t is the cosmic look-back time, H 0 is the Hubble 
constant, and rj is a variable scaling parameter. 

3. Quasar-like density evolution, where the burst rate 
density is proportional to exp [ — (z — z peak ) 2 /2<r 2 ], i.e., 
centered on a redshift z peak with a spread of a. 

Of course one has an infinite number of a priori choices for 
n c (z), however we choose these three for their simplicity so that 
the effects of PDE can be more easily understood in the 
absence of many free parameters. 

In the following section, we describe the analysis technique 
used for each of these three evolutionary scenarios. Subsequent 
sections of the paper discuss results from each of these three 
scenarios in detail. In the final section of the paper, we discuss 
conclusions regarding the modes of PDE compatible with the 
BATSE data. 


2. ANALYSIS TECHNIQUE 

The methodology by which compatibility of a given evolu- 
tionary scenario is determined in this work is identical to the 
analysis of Emslie & Horack (1994). For a cosmological burst 
distribution, the normalized integral moments {N im) } of the 
differential BATSE brightness distribution N(P) can be 
expressed as (see Horack & Emslie 1994) 


W-H = fo"“ n c (z)f(z)r(z) 2 - 2 "(l + z)-™~ l dz 

1 ' 1 \r"n c (z){\ + zy'f(z)r{zydz ’ 


( 1 ) 


where 


J [>o(l + z) 3 + (1 + qo ~ 3<7 0 )(1 + z) 2 + (<x 0 - <?o)] 1/2 ' 

( 2 ) 

Here, <r 0 is the density parameter, q 0 is the deceleration 
parameter, and a is the burst power-law photon spectral index, 
which we take to be a representative value of 2 throughout this 
analysis. The {0£°} are the normalized integral moments of the 
function <£*(L) which represents the number of bursts with 
peak luminosities between L and L + dL that are accessible to 
the BATSE detector, and also incorporates various normal- 
izations such as the current value of the scale-factor S 0 and the 
Hubble constant H 0 . Although the relationship between 
and <p(L), the luminosity distribution of all bursts, depends on 


the form of n c (z), for the purposes of this analysis, one only 
needs to assert that the function <t> m (L) exists and that it is 
non-negative everywhere. 

The integral moments of the BATSE differential N(P) dis- 
tribution, along with their associated uncertainties, are easily 
measured from the observational data and are represented on 
the left-hand side of equation (1). Horack & Emslie (1994) 
provide a thorough discussion of how these numbers are 
obtained. Briefly, the mth integral moment {N (m) } is simply 
computed as 

{N (m) } = <P m > , (3) 

where P is the 0.256 ms peak flux of gamma-ray bursts in the 
energy range 50-300 keV. Only those bursts with peak fluxes 
in excess of 0.5 photons cm -2 s -1 are utilized, to avoid the 
instrumental selection effects that are known to bias the inten- 
sity distribution below this value (Fishman et al. 1994). Uncer- 
tainties in these quantities, like those in the computation of the 
angular distribution moments (cos 0} and (sin 2 b>, are due to 
both uncertainties in the peak flux P of each individual 
gamma-ray burst, as well as statistical uncertainties due to the 
finite number of bursts in the sample. Horack et al. (1994a), 
Briggs et al. (1994), and others have shown that in computing 
the angular distribution moments, the total uncertainty in the 
measured value is dominated by the 1 / v /n statistical uncer- 
tainties due to the finite number of bursts. Not surprisingly, 
this is also the case with the brightness distribution moments 
{AT , ”’ > } (Horack & Emslie 1994; Pendleton et al. 1995). 

After computation of the left-hand side of equation (1) from 
the data, if one subsequently assumes a cosmological model 
(tr 0 , q 0 ) and a burst distribution n c (z) out to a limiting redshift 
z mak , the integral on the right-hand side of equation (1) can be 
computed and the set of numbers { <l>! k m) } and their uncertainties 
can be derived. If the assumed cosmological scenario is valid, 
these numbers are the integral moments of the associated 
</>„(L) function, from which minimal information on <t>„{L) and, 
since n c (z) has been assumed for the given model, <p(L) can be 
deduced. Whether or not this information is accurate depends 
on the validity of the assumed cosmological scenario. 

It might seem that, with no a priori information on either 
4>(L) or n c (z), little information can in fact be obtained through 
such methods. However, the integral moments of a given func- 
tion are not all independent quantities. For example, for a 
function /(r), the second moment <r 2 > must be larger than (or 
equal to) the square of the first moment <r> 2 , in order that/(r) 
have a realizable (non-negative) variance. Consequently, an 
arbitrary set of computed numbers need not be valid moments 
for the {(p 1 ™’}. In order that a given set of numbers represent 
the moments of a positive semi-definite function, these 
numbers must satisfy conditions known as the Wald inequal- 
ities (Wald 1939). These inequalities state that a set of numbers 
M r , Ms> Mr can be realized as moments of order r,s,t(r < s < t) of 
some non-negative distribution function defined on the inter- 
val [0, oo) if and only if 

Mr>0; H S >H S J T \ Mr Ms Mr 2: M?mX • (4) 

For the problem at hand, then, if a set of numbers {<t> , ” ) } are 
obtained through the assumption of a cosmological model and 
employment of the BATSE N(P) normalized integral moments 
{N (m) }, and these derived numbers violate any of the inequal- 
ities in equation (4), not only can these numbers not be the 
moments of the $*(L) function, they cannot be the moments of 
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any function. One is therefore forced to conclude that the 
cosmological model and evolutionary scenario that was 
chosen to generate these numbers are incompatible with the 
BATSE data for all ip(L) at the appropriate level of significance 
indicated by the propagated uncertainties. In this manner, the 
viability of certain cosmological evolutionary models can be 
ruled out. 

Specifically for this analysis, we have chosen those moments 
from the observed N(P) distribution which can be measured to 
an accuracy of ~ 10%— 15% or better, corresponding to values 
of m = — 1, — 2, ... , — 6. We then select values in the (<r 0 , q 0 ) 
plane, combined with an assumed form for n c (z), and use the 
wald inequalities to determine the largest limiting redshift z max 
for which a mathematically tenable function can be con- 
structed. This generates a volume of parameter space in which 
all viable cosmological scenarios with the assumed comoving 
density function must exist. Figure 1 shows this volume as 
obtained by the method of Emslie & Horack (1994) for a non- 
evolving comoving density n c (z) = constant. Emslie and 
Horack also have computed surfaces corresponding to the + 1 
a and +3 a uncertainties applied to the observed moments 
The basic form of these surfaces is the same. 

In addition to constraining the possible cosmological 
models consistent with the BATSE data, further information 
may be obtained by realizing that the upper surface to this 
parameter volume consists of cosmological models which 
admit luminosity functions <t> t (L) which are just compatible 
with the Wald inequalities, i.e., »<5(L — L 0 ). The height of the 
upper surface therefore corresponds to the maximum redshift 
obtained from the BATSE data for a given [n c (z), <t 0i <? 0 ] 
assumption with a burst population that is mono-luminous. 
The farther a model lies below this limiting surface, the greater 
the flexibility in choosing a 0*(L). 

For example, if z max <? 1, the nonevolving bursts are embed- 
ded in a Euclidean space and evolutionary effects are negligi- 
ble. The deviation from a — 3/2 power law in the log j V( > P ) 
profile is then due to the instrumental truncation of low- 
luminosity bursts at large distances (which therefore have a 
brightness below the BATSE threshold), and the rate at which 
the curve deviates from the — 3/2 power law is a reflection of 
the shape of the luminosity distribution. Steep 0 # (L ) cause a 
relatively rapid fall-off as is observed, while shallow 0*(L) 
would cause a more gradual fall-off in contradiction to the 



Fig. 1. — Limiting redshift z„„((t 0 , <? 0 ) derived for a nonevolving population 
of cosmological gamma-ray bursts. Cosmological scenarios that are compat- 
ible with the BATSE data lie on or below the surface. (Emslie & Horack 1994). 


BATSE observations. The restriction on overly shallow forms 
of <f>+(L) is consistent with the low height of the z max surface, 
which does not permit models to lie substantially below it. 

Therefore, regions of the parameter volume which prohibit 
either large values of z max or which prohibit the model from 
lying substantially below the upper surface are therefore con- 
strained to have relatively steep and/or narrow observed lumi- 
nosity ranges (see Horack, Emslie, & Meegan 1994b; Emslie & 
Horack 1994), depending on the form of the radial distribution. 
As can be seen from Figure 1, this is indeed the case for most 
nonevolving scenarios with q 0 larger than zero, including all 
models for which A = 0. We shall see in subsequent sections 
that evolutionary models can permit a larger value of z max and 
hence a shallower form of <£*(L). 

3. POWER-LAW EVOLUTION IN REDSHIFT 

The first evolutionary scenario we examine consists of a 
comoving burst density that evolves as 

n c (z) cc(l + zf , (5) 

where P is a variable index that can be used to increase or 
decrease the amount and direction of the evolution present in 
the scenario. Although not necessarily motivated from observ- 
ational data, this form of evolution is useful in determining 
what effects are to be found for various amounts of evolution 
in the density of bursts. 

We have determined the largest permissible value of 
z max(°o> 9o) f°r a variety of values P ranging from — 2 to 2. The 
<7o), Zmax + 1 and z max + 3 a surfaces with p = 0 are 
obviously identical to those of the nonevolving case (§ 2, Fig. 1, 
Emslie & Horack 1994). For comparison with a moderately 
evolving n c (z), Figure 2 shows the derived z max (<r 0 , q 0 ) limiting 
surface for p - 1.7; corresponding to nearly 3.25 times the rate 
of bursts per unit volume of comoving space at z = 1 than are 
present locally. The general shape of the two limiting surfaces 
in Figure 1 and Figure 2 are the same, however with the 
increased value of p, larger values of z max are admitted and the 
plateau region adopts a slightly more curved shape. 

This result is easily understood when one considers the log 
J^( > P) versus log P curve. With an excess of bursts at large 
redshifts/distances, the curve at the low-P end of the distribu- 
tion tends to display a steeper slope than would be found for a 
population with no evolution. In fact, in a Euclidean distribu- 
tion, the slope would be steeper than the uniform density value 



Fig. 1 — Limiting redshift q 0 ) for bursts evolving as n ( (z) cc (1 

+ z)'\ 
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of —3/2. Larger redshifts are therefore required over the 
observed range of intensities to introduce stronger cosmo- 
logical bending of the jV(>P) curve to the flatter (slope 
~ —0.8) form observed in the BATSE data. This requirement 
of larger z max is therefore a general one for all evolutionary 
scenarios which have monotonically increasing n c (z) and a 
narrow 4> t (L). 

Figure 3 shows plots of the maximum limiting redshift z^ x 
determined for three simple cosmological scenarios as a func- 
tion of the index /?. All three of these scenarios have a zero 
cosmological constant, with values of the deceleration param- 
eter q 0 equal to 0 (diamonds), 0.5 (triangles), and 1.0 (squares). 
Cosmological models with A = 0 and power-law density evol- 
ution must possess limiting redshifts equal to or less than the 
value of z max shown in order to be consistent with the observed 
BATSE brightness distribution. Furthermore, models with 
limiting redshifts very near the value of z max in the figure corre- 
spond to nearly mono-luminous bursts. 

This figure permits one to draw various conclusions regard- 
ing the amount of power-law source evolution that can be 
present in the burst population depending on the actual limit- 
ing redshift in the BATSE burst sample, and some general 
assumptions regarding the type of luminosity function the 
observed bursts possess. 

As an example, one independent measurement of the limit- 
ing redshift of bursts in the BATSE sample may have been 
given by the work of Norris et al. (1995), who used the duration 
distributions of bright and dim bursts to search for cosmo- 
logical time-dilation. It is well-known that the average dura- 
tion for a set of gamma-ray bursts at redshift z dim , when 
compared with the average duration from a set of bursts at 
smaller redshift z br „ can be expressed as 

<T iim > _ 1 + ^dim 

<T brt > ” 1 + z brt ' U 

Norris et al. have measured this ratio to be slightly in excess of 
2. Combined with the assumption that the bright gamma-ray 
bursts are located at a redshift of z brt « 0.3 (Norris et al. 1995), 
their analysis indicates that the BATSE data display time- 
dilation effects consistent with the dimmest bursts at limiting 
redshifts of ~2. A similar time-dilation effect has also been 
observed by Wijers & Paczynski (1994). However, Emslie & 
Horack (1994) have shown that a limiting redshift of z max % 2 is 



Power— Low Evolution Scole Factor (3 

Fig. 3. — Limiting redshift vs. power-law evolution index /} for three 
A = 0 cosmologies with an evolving burst population n c (z) oc (1 + zf. 


inconsistent with the BATSE data and a nonevolving popu- 
lation of gamma-ray bursts at the +3 a level. Evolution in the 
source population provides a means of reconciling the appar- 
ent limiting redshift obtained with time-dilation measurements 
and the BATSE N(P) distribution. 

If one chooses this value of z lim — ^ as an independent mea- 
surement of the actual limiting redshift in the BATSE data, 
evolutionary scenarios that are not inconsistent with the 
BATSE data at some level can only be found in Figure 3 where 
the model lies near z max = 2, and below the curves displayed in 
the figure. This restricts the admissible forms of power-law 
evolution functions to those where ff is larger than 1. 4- 1.6, 
depending only slightly on the chosen value of q 0 . Therefore, if 
the Norris et al. (1994) limiting redshift of 2 is the correct value 
for the BATSE data, power-law evolutionary scenarios consis- 
tent with the BATSE N(P) distribution and the measured limit- 
ing redshift must be at least as strongly increasing as n c (z) oc (1 
4- z) 15 , and negative values of / ? are excluded. This conclusion 
is furthermore independent of the luminosity function chosen 
for the bursts. 

An additional assumption regarding the range of luminosity 
present in the BATSE data can constrain the maximum 
amount of evolution that is permissible. For example, if one 
additionally assumes that the luminosity function of the 
observed gamma-ray bursts is rather narrow (an assumption 
that may be supported in the context of A = 0 cosmologies by 
the rapid turnover in the integral number-intensity distribu- 
tion observed by BATSE; Emslie & Horack 1994; Mao & 
Paczynski 1993), one further constrains the locus of points in 
Figure 3 consistent both with the BATSE N(P) distribution 
and these additional studies. Consistent points now not only 
must lie very near to the plotted curves in Figure 3. For 
an absolute redshift limit of z ^ ~ 2 and standard candle 
burst sources, the amount of power-law PDE in the co- 
moving density cannot be any more or less severe than about 
P w 1. 5-2.0. 

4. EXPONENTIAL EVOLUTION IN COSMIC LOOK-BACK TIME 

The second form of PDE is exponential, with the comoving 
density of bursts following the form 

« c (z) oc exp [r/r(z)H 0 ] , (7) 

where t (z) is the cosmic look-back time given by 


=-T 

H 0 Jo 


(1 + 2') 


with H 0 being the Hubble constant, /(z) is given by equation 
(2), and q is an exponential scaling parameter. 

The procedure for analysis of this model of evolution is 
formally identical to the previous section. The z max (<T 0 , 9o) 
surface is constructed using the moments of the BATSE N(P) 
distribution along with the Wald inequalities. For points 
above the surface, no physical <f>+(L) function can be con- 
structed given the BATSE data and the assumed cosmological 
parameters. Figure 4 shows two of the z mzx (a 0 , q 0 ) surfaces for 
values of q equal to 3 and 9. 

Qualitatively, these surfaces are similar in shape to the limit- 
ing surface derived with no density evolution (Fig. 1). As was 
seen in the case of power-law density evolution, as the amount 
of evolution in the burst population increases (q becomes 
larger), the maximum limiting redshift also increases to provide 
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Fig. 4. — Limiting redshift z„„(<r 0 , q 0 ) surfaces for bursts evolving as tt,(z) oc exp (tyiH 0 ): (a) t] = 3; (ft) rj = 9 


additional flattening of the log Jf(> P) distribution that offsets 
the increased number of bursts at larger distances. 

In Figure 5, we present the maximum allowable redshift as a 
function of the exponential scale factor rj for three cosmo- 
logical models with A = 0 and q 0 =0, 0.5, and 1.0. Again, 
cosmological scenarios with A = 0 and exponential density 
evolution that are consistent with the BATSE observation of 
N(P) must have limiting redshifts with values less than or equal 
to the values of z max shown in the figure. As before, the plotted 
points correspond to consistent scenarios with standard- 
candle burst populations. 

Qualitatively, the curves shown in Figure 5 for exponential 
source evolution are similar to those shown previously for the 
power-law evolutionary mode, and similar conclusions to § 3 
may be drawn from Figure 5. 

5. QUASAR-LIKE GAUSSIAN DENSITY EVOLUTION 

The two previous forms of PDE explored in this paper are 
both rather ad hoc, with little or no physical underpinning to 
rationalize their use as possible forms of source evolution for 
gamma-ray bursts. However, as mathematical functions, they 
offer insight into the types of effects that source evolution can 
have on the observational data. Both of the previous functions 
are monotonically increasing with redshift and/or look-back 


time. While it is possible that over the range of observable 
redshifts the comoving burst rate density does continue to 
increase, it is reasonable to assert that at some redshift or 
look-back time, this comoving density of bursts may decrease 
to zero, certainly by the time that one reaches the bigbang. It is 
also possible that, as in the case with quasars, the density of 
sources does begin to decrease at a distance/redshift that is 
interior to the maximum being sampled by the instrument. A 
source density with a turnover interior to the volume of space 
sampled by the experiment will have very different conse- 
quences for the observed N(P) distribution, compared to a 
source density that continues to increase out to the experi- 
ment’s threshold. 

Quasars are among the earliest visible objects to form in the 
universe, currently detectable to redshifts of order 5 (Schneider 
et al. 1989). As such, they may be the first indicators of the 
formation of large-scale structures, with evolutionary connec- 
tions to AGNs, BL Lac objects, Seyfert galaxies, and “ normal ” 
galaxies that we observe in nearby space. Figure 6 is adapted 
from Hartwick & Schade (1991) and contains the integrated 
space density of luminous quasars (absolute magnitude M B < 
— 26) for the two parameter sets q 0 = 0.1, n q = — 1 ( open 
circles) and q 0 = 0.5, a, = —0.5 (filled circles ), where a, is the 
power-law spectral index of the quasars. 



Fig. 5. — Limiting redshift z vs. exponential scale-factor tj for three A = 0 

cosmologies. 



Fig. 6. — Measured comoving number density of quasars for two different 
assumed spectral indices and cosmologies. Adapted from Hartwick & Schade 
(1991). 
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Although difficult to determine precisely, the comoving 
density of quasars increases with redshift out to a broad 
maximum around z « 2-2.5, followed by a general decrease in 
the comoving density of sources. The solid curve shown in the 
figure is a simple Gaussian representation of the filled-circle 
data, with a mathematical form 

Hc (z)oeexp[ ~ (2 ~ Z / cak) -] (9) 

with values of z peak = 2.15 and a = 0.45. This simple curve 
does not fit the data precisely, however it does capture the 
general characteristics of the source density increase, followed 
by a broad maximum and then a decline. Additionally, this 
functional form can also yield a reasonable representation of 
the open-circle data (q 0 = 0.1, a, = -1) simply by adjusting 
the amplitude and leaving the values of z peak and a fixed. 

Gamma-ray bursts at cosmological distances require the 
existence of some type of progenitor object which itself may 
follow a similar evolution in comoving density to that of the 
quasars. This notion may be especially relevant for objects like 
AGNs, which are observed to display an evolutionary connec- 
tion to the quasars, and have been speculated to be one of the 
possible sites of gamma-ray burst production (Brainerd 1994). 
The redshift at which the peak density occurs, however, may be 
different for the gamma-ray bursts than shown in Figure 6 
because of a “ lag ” in cosmic time between the formation of the 
quasars and the development of the burst progenitor objects. 
Similar lags occur in NS-NS merger scenarios, where binary 
formation and the merger event would have a large range of 
lag times, depending on the initial binary separation. 

We therefore have chosen to explore the effects of PDE in 
the context of the BATSE gamma-ray burst data utilizing the 
comoving number density of Equation (9) with various values 
of z peak and a. The motivation to use such a function comes 
from both its mathematical nature, with an actual decrease in 
the comoving density of sources beyond a certain redshift, as 
well as for its possible physical implications regarding the evol- 
ution of objects known to exist at cosmological distances. 

Figure 7 displays the z max (<7 0 , q 0 ) surface obtained with a 
comoving number density prescribed by equation (9), with a 
value of z peak =1.0 and a width a = 0.45. This figure is again 
qualitatively similar to the previous surfaces shown with other 
forms of density evolution. Despite the similarities, there are 
some interesting differences as well. The plateau region for this 
type of evolution appears flatter, with even less z max depen- 



Fig. 7. — Limiting redshift q 0 ) computed for a burst population 

which evolves according to eq. (7). For this figure, = 1.0 and a = 0.45. 


dence on the exact values of a 0 and q 0 for these cosmologies. 
For smaller values of z peak , the plateau region is more extensive 
in area and requires smaller values of z^ as z peak is decreased. 
For z peak equal to zero (i.e., the maximum comoving density of 
bursts occurs locally), a maximum limiting redshift of z max « 
0.5 is found for all cosmologies with q 0 > — 1, regardless of the 
value of <7 0 . This is because the fall-off in n c (z) is already suffi- 
cient to flatten the — 3/2 Euclidean power law to the observed 
form without the need for significant cosmological effects. At 
large values of z peak > 2, there is no obtainable constraint on 
z mal , as most of the bursts are concentrated in a relatively 
narrow range of redshift (a “shell-model”), and an assumed 
luminosity function similar in form to the observed brightness 
distribution can always be found to make the observed spatial 
distribution produce a brightness distribution that is consis- 
tent with what is actually being observed by BATSE. 

Figure 8 contains plots of z max versus z peak for a = 0.45 and 
cosmologies with q 0 = 0, 0.5, and 1.0, each with a zero cosmo- 
logical constant. Both the measured value of z max is shown 
( open symbols) as well as the + 3 a upper limit on z max ( filled 
symbols). As with the previous two similar figures (Figs. 3 and 
5), brightness distributions consistent with the BATSE data 
can be found for a given z peak and the PDE shown in equation 
(9) provided that a suitable luminosity function is chosen and 
the limiting redshift is less than or equal to z max . Points plotted 
in the figure correspond to models with nearly mono-luminous 
burst populations, and the required luminosity function 
becomes broader as the limiting redshift is decreased below the 
plotted value of z max . 

Interesting conclusions can be drawn from Figure 8 with an 
additional assumption regarding the actual limiting redshift in 
the BATSE data. If we again choose the value of z lim « 2, it is 
clear from Figure 8 that the value of z peak cannot be less than 
~ 1.4 with a 3 a lower limit of ~ 1.0, regardless of the value of 
q 0 chosen. Given the measurement of Norris et al. (1994), it is 
impossible for a quasar-like comoving burst density of the 
form of equation (9) to peak at a redshift smaller than z peak ~ 
1.0 and still remain consistent with the observed BATSE data 
in a A = 0 universe. 

The effect of modifying the width of the gaussian distribu- 
tion ( <t ) in equation (9) is shown in Figure 9 where the 
maximum limiting redshift z max is plotted as a function of 1 /ct. 
For large values of z peak , making the distribution narrower 
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Fig. 9. — Limiting redshift z rnxx as a function of the inverse of the Gaussian 
njz) width a. The computation is artificially truncated at z mxx = 10. 


(increasing 1 /<r) increases the maximum limiting redshift. The 
narrower distribution more closely approximates the “shell 
model” discussed earlier, and is therefore compliant with the 
data over a larger range of limiting redshifts. If z peak is small, 
however, the limiting redshift is decreased as the distribution is 
made narrower. In this case, the steep drop-off in the burst rate 
density is itself sufficient to cause the observed bending in the 
log Jf( > P) distribution without the introduction of significant 
cosmological effects that would be introduced by observing to 
large redshifts. 

As the distribution of burst rate density is given a larger 
width (making l/a smaller), it more closely approximates the 
case of a distribution that is roughly constant over a large 
range of redshifts, a scenario that has been widely studied in 
the literature (Wickramasinghe et al. 1993; Mao & Paczynski 
1992; Emslie & Horack 1994). In agreement with these pre- 
vious studies, the limiting redshifts in these cases approach 
z max ~ 1 as a becomes large. 

If we again adopt z max = z Um = 2, as suggested by the work 
of Norris et al. (1995), then for monoluminous bursts we see 
that z peak must be greater than about unity, and that the width 
of the gaussian distribution must be at least as wide as a « 0.3 
depending on the exact value of z peak . This is a rather pleasing 
result; such an evolutionary scenario creates an enhancement 
of bursts at just the right redshift, and over just about the right 
range of redshifts to produce the quick roll-over observed in 
the BATSE log > P) distribution. This rapid roll-over has 
been interpreted, in a Euclidean geometry, by a width in the 
observed burst luminosity distribution of about a factor of 6 
(Horack et al. 1994b). For z pcak = 1 and a — 0.33, the spread in 
redshifts within which 80% of the bursts are contained is from 
0.58 to 1.42, or a range of a factor of 2.45. The square of this 
quantity is indeed about 6, the exact range of L-values within 
which 80% of the bursts must be contained if a Euclidean 
scenario is valid (Horack et al. 1994b). 

6. DISCUSSION 

The technique of integral moments has provided insight into 
three evolutionary scenarios for the comoving rate number 
density of gamma-ray bursts. In all three of the scenarios that 
have been examined, power-law evolution, exponential evolu- 
tion, and a Gaussian quasar-like evolution function, many of 
the results are qualitatively similar. We observe for a large 


number of (<x 0 , q 0 ), including all A = 0 models, that there exists 
a broad z m .Jo n , q 0 ) plateau that limits the maximum possible 
observed redshift. As the amount of evolution is increased in a 
given model, the admissible value of z max increases. 

Introducing an externally derived value for the limiting red- 
shift in the BATSE sample provides a means by which the 
amount of evolution in the burst population can be con- 
strained. By employing the limiting redshift Zlim ~ 2 found by 
Norris et al. (1995) through the analysis of temporal burst 
profiles from the BATSE data, all mathematically tenable evo- 
lutionary scenarios with limiting redshifts that are smaller than 
this value for all possible luminosity functions are excluded. 

A special case is the scenario in which the comoving burst 
density is a constant, i.e., no evolution. (See Figs. 1 and 3 with 
/? = 0, Fig. 5 with t] = 0, or Fig. 9 with 1/tr = 0.) Very few of the 
nonevolving cosmological models, indeed only accelerating 
models where q 0 < 0, can have limiting redshifts of 2, even if 
one chooses to use the + 3 a upper limits to the derived z max - 
value as the quantity of interest. On the basis of this analysis 
and the strict application of an independent measurement of 
the limiting redshift in the BATSE data, we therefore conclude 
that if the actual limiting redshift in the BATSE data is z lim ~ 
2, then there must be some form of evolution present in the 
gamma-ray burst population unless A > 0. Furthermore, this 
evolution must be of a form to either increase the comoving 
density with redshift, increase the luminosity of bursts with 
redshift, or some combination of both. 

If the evolution is PDE, for A = 0 cosmologies, power-law 
evolution must display a rate density increase at least as fast as 
(1 + z) 1 5 , and exponential evolution must increase at least as 
fast as exp (3 tH 0 ) to avoid inconsistency with the BATSE N(P) 
distribution at some level. A gaussian-like quasar evolution 
function with a a ~ 0.45 must peak at a redshift no less than 
z ~ 1.0. Employing different values of q 0 in these A = 0 
cosmologies makes only a minor difference in these results (see 
Rutledge, Hui, & Lewin 1995). 

The necessity of evolution in cosmological scenarios has 
been discussed previously in the literature by Fenimore et al. 
(1992) who utilized early BATSE data in conjunction with data 
from PVO. However, the form of evolution that was stated to 
be required in their analysis was in the opposite sense to that 
found here, with fewer or fainter bursts at large redshifts. A 
re-examination of the effective PVO lifetime and subsequent 
analysis of the combined data (Fenimore et al. 1993) retracted 
the previous requirement of evolution, showing that the 
BATSE-PVO brightness distribution is consistent after all with 
a nonevolving population to limiting redshifts of order unity. 

The limiting z raax (cr 0 , q 0 ) surface derived for the three evolu- 
tionary scenarios corresponds to a mono-luminous population 
of gamma-ray bursts. Therefore, if one further requires that 
bursts possess identical luminosities in addition to the previous 
constraint that z max ~ 2 in the BATSE data, one severely 
restricts the amount of PDE source evolution that can be 
present in the data. Again for A = 0 cosmologies, power-law 
type evolution of the form n c (z ) oc (1 + zY is incompatible with 
the observational data and these assumptions outside the 
range 1.5 < P < 2.0, depending somewhat on the value of q 0 
chosen. In the case of exponential evolution with cosmic look- 
back time, the comoving density function n c (z) oc exp (rjr H 0 ) is 
confined to the region 3 < tj < 6. 

With an assumed mono-luminous burst population and an 
actual limiting redshift of 2, the quasar-like Gaussian evolution 
function of Hartwick & Schade (1991) shown in Figure 6, must 
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peak somewhere between 1.0 < z peak < 1.5. Evolutionary 
density functions of this type which peak at redshifts smaller 
than z = 1 are incompatible with the BATSE data regardless 
of one’s assumption of the luminosity function, provided that a 
limiting redshift of z = 2 is present in the BATSE N(P) dis- 
tribution. 

In this work, we have focused only on PDE modes of evolu- 
tion because of their mathematical simplicity, their possible 
relation to certain classes of models, and the applicability of 
integral moment analyses to these forms of evolving burst 
populations. There remains the possibility that the luminosity 
function of the bursts also evolves. Although the integral 
moment analysis technique is not equipped to handle insepara- 
ble forms of <t>(L , z), a qualitative understanding of luminosity 
evolution and its effects on the observed N(P) distribution can 
be understood on the basis of this analysis. Luminosity evolu- 
tion which provides for more luminous bursts at large redshifts 
with no concurrent source evolution will require the presence 
of larger limiting redshifts in order that the resulting ^ V(>P ) 
distribution be consistent with the observed data. Alternative- 
ly, a decrease in burst luminosity with increasing redshift will 
work in the opposite direction, reducing the allowable limiting 
redshifts consistent with the BATSE data. 
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Given the freedom to invoke both luminosity and source 
density evolution, an extremely large number of combinations 
will be able to generate a theoretical brightness distribution 
that mimics the BATSE observed N(P) distribution. However 
in many cases, especially those scenarios with only moderate 
limiting redshifts, these two evolutionary modes must work 
against each other to produce a resulting intensity distribution 
consistent with the BATSE data. Cosmological models with 
very strong evolution in the luminosity function that is 
compensated by strong evolution in the source density func- 
tion (in the opposite sense) in order to produce a brightness 
distribution that is consistent with the BATSE data, itself 
consistent with a nonevolving, mono-luminous population 
(Wickramasinghe et al. 1993; Fenimore et al. 1993; Emslie & 
Horack 1994) would be suspect through an Occam’s Razor 
argument. 
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